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ABSTRACT 

The fragmentation of shocked flows in a thermaUy bistable medium provides a nat- 
ural mechanism to form turbulent cold clouds as precursors to molecular clouds. Yet 
because of the large density and temperature differences and the range of dynami- 
cal scales involved, following this process with numerical simulations is challenging. 
We compare two-dimensional simulations of flow-driven cloud formation without self- 
gravity, using the Lagrangian Smoothed Particle Hydrodynamics (SPH) code VINE 
and the Eulerian grid code Proteus. Results are qualitatively similar for both methods, 
yet the variable spatial resolution of the SPH method leads to smaller fragments and 
thinner filaments, rendering the overall morphologies different. Thermal and hydro- 
dynamical instabilities lead to rapid cooling and fragmentation into cold clumps with 
temperatures below 300 K. For clumps more massive than 1 Mq pc~^, the clump 
mass function has an average slope of —0.8. The internal velocity dispersion of the 
clumps is nearly an order of magnitude smaller than their relative motion, rendering 
it subsonic with respect to the internal sound speed of the clumps, but supersonic as 
seen by an external observer. For the SPH simulations most of the cold gas resides at 
temperatures below 100 K, while the grid-based models show an additional, substan- 
tial component between 100 and 300 K. Independently of the numerical method our 
models confirm that converging flows of warm neutral gas fragment rapidly and form 
high-density, low-temperature clumps as possible seeds for star formation. 
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1 INTRODUCTION 

The highly filamentary morphology of molecular clouds 

I MCs) and their observ e d non-therma l line -widths 
Falgarone fc PhiUipsI Il990l : IWilliams et~al] |2000') point 
to MCs being highly dynamical objects. Observational 
evidence suggests that star formation in local MCs such 
as Taurus is rapid once molecular gas is available, and 
that the parental clouds are short-lived (jElmegreenl I2OO0I : 
iHartmann et all I2OOII : iHartmannl |2003| ). The dynamical, 
or turbulent, nature of MCs is assumed to play a crucial 
role in the process of s tar formation via turbule nt frag- 
mentation l|Larsonlll98ll : [Mac Low fc Klessenll2004l . see also 
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iBallesteros- Paredeslliooel . for a summary of the effects and 
interpretation of turbulence in MCs). 



Because of the rapid onset of star formation, the cloud 
formation process needs to provide the MC with the ob- 
served turbulence and substructure. Moreover, global geom- 
etry and gravity considerations mandate that this substruc- 
ture be non-linear (Burkcrt & Hartmann 20041 ). i.e. a phys- 
ical process is needed that can imprint non-linear density 
perturbations in the proto-cloud during its formation. These 
requirements have led to the scenario of flow-driven cloud 
formation, where molecular clouds are assembled by large- 
scale converging flows of atomic hydro gen corresponding to 
the warm neutral medium |B allcstcro s-Paredes et al.l 19991: 
IHartman n et al. l l200ll : IHartmann.. 20031 . see also lElmegreenI 
1993, 200(1). The rapid fragmentation is driven by a combi- 
nation of strong thermal and dynamical instabilities, dom- 
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inated by the therm al instabihty (TI, iFieldl [l965l . see also 
iHeitsch et al.ll200^ for a discussion of timescales) . 

Large-scale gas flows are ubiquitous in the Galaxy. 
They might be driven loca ll y by supernova explo- 
sions (lElmegreen fc Ladal 1 19771 : iMcCrav fc Kafatoj 1 19871 : 
iNigra et al.ll2008l ) or globally by shear motions i n the Galac- 
tic d isk, by global graviation al instabilities jYang et al. | 
2007[) . gas infall from the halo (|Lacey fcF all 1985: Mirabel 
1983), interactions with the central bar (jRobcrts ct al..,197S : 



Combes fc Gerinlll985l ) or satellite galaxies. On extragalac- 
tic scales collisi ons of galactic dis ks trigger gas flows, shocks 
and starbursts (|Mihos fc Herngu ist 1996; Naab ot al. 2006; 
iKarl et al.|[201ol ). 

High-resolution simulations in two 



ll Audit fc Hennebelld |2005| ; [Hennebelle et all l2007l : 



iHennebeUe fc AuditI l2007l: iHe itsch et al.' '2005' 
and three dimensio n s Ivazauez-Semadeni et al 



2006") 
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20091 ) have demonstrated that the flow-driven formation 



of clouds is indeed a natural and elegant way to provide 
the clouds with the observed turbul e nce a nd substructure. 
Except for IVazguez-Semadeni et al] (|2007l ). the above au- 
thors used various grid-based methods. However, due to the 
high density contrast and the strong TI in the intermediate 
temperature regime between 300 < T < 5000 K, the spatial 
scales of the forming cold filaments and clumps shrink 
dramatically. This problem becomes even more severe with 
gravity acting upon the cold regions. 

Due to its Lagrangian nature, the strength of Smoothed 
Particle Hydrodynamics (SPH, e.g. Monaghan 1992) re- 
sides in its capability to follow fluid flows at "arbi- 
trary" spatial resolution. In particular, dissipative prop- 
erties of SPH methods do not depend on geometry or 
direction, which can cause issues for grid-based met hod s 
without physical dissipation control ((Rampp et al.i 1 19981 ). 
SPH has been used to model the inte rstellar medium o n 
all scales, from planet formation (e.g. iMaver et al.ll2003 ). 
formation o f protostella r disk s from rotating and turbu- 
lent cores fWalch et al.' '2009", '2010'), formation of stars 
and cores (IKlessen .1997 ; Bate et al. 2003), ionizatio n feed- 
back from massive stars (jGritschneder et all |2009| . I2OI0I : 
iBisbas et all |2009|) formation of MCs in g alactic disks 
dDobbs et al.l I2OO6I : iDobbs fc Bonnelll |2007l: iDobbs et al ' 
' 20081' 



;alaxy evolution (ISteinmctz 1996; Hernguist fc Katz 



19891 : l^ringel 2000l : iNaab et al. 2006. 200"7l). to large scale 



simu la tion of galaxy formation ( e.g. Springel fc Herngui"stl 
I2OO3I ). IVazguez-Semadeni et aO 1)2007^ present models of 
flow-driven MC formation using SPH. However, SPH has its 
own inherent limitations, among others a limited mass res- 
olution, and spatially varying dissipative properties. Hence 
the question remains whether SPH and grid-based methods 
give similar results for a specific astrophysical problem, and 
what the method-intrinsic differences are. iNagamine et al.l 
l|2005l) showed for cosmological applications that both ap- 
roach es give statistically similar results, and lKlessen et ahl 
2OOOI ) concluded the same for models of self-gravitating 
driven turbulence in a period ic box of isothermal gas, al- 
though Kitsiona s et al.l (|2009| ) show for the latter applica- 
tion that SPH is more dissipative than grid models for de- 
caying turbulence. 

Yet it remains unclear whether this agreement between 
methods holds for models of a thermally unstable, turbu- 



lent interstel l ar m ediu m. This is of part icular interest as 
lAgertz et all (|2007l ) and lJunk et all l|2010l ) have pointed out 
SPH's limitations evolving shear flow instabilities (see, how- 
ever, [Price 2008; R ead et al. 2010; Abo] 2010) - a crucial 
ingredient in the formation o f turb ulent clouds from collid- 
ing flows (|Heitsch et al.ll2006l . [200^ ). 

This paper compares the Eulerian and Lagrangian ap- 
proach to the problem of flow-driven cloud formation in 
two dimensions. The co ncept of the simulations follows 
IHeitsch et al] (|2005l . I2OO6I ) and is summarized in g2] We find 
(© ^ good agreement between SPH and grid-based results 
in terms of the rapid onset of fragmentation and, to a lesser 
degree, of morphology. 

The total amount of mass accumulated in the cold gas 
component (T < 300 K) is smaller in the SPH models, yet 
the temperature distribution within the cold gas is concen- 
trated towards lower temperatures. The grid models show 
a wider spread of temperatures, but an overall higher mass 
in cold gas. We attribute this di fference to the over- cooling 
problem in SPH described by e.g. lPearce et al.l (|l999l ). Line- 
of-sight velocity dispersions as seen by an observer are gener- 
ally supersonic with respect to the cold gas, yet, the disper- 
sions within the cores a re approximately s onic, in agreement 
with observations (e.g. I Andre et al.l 120071 ). The strength of 
the Lagrangian approach is demonstrated when deriving the 
mass distribution of cold dense fragments forming in the col- 
liding flows. 



2 THE MODELS 

All SPH simulations ar e run with the para l lel shared mem- 
ory S PH code VINE l|Nelson et all l2009l : IWetzstein erall 
[200y). The code solves the asymmetric form of the energy 
equation and uses the standar d implementatio n of the arti- 
ficial viscosity to treat shocks (|Monaghanlll992l ). SPH parti- 
cles have individual time-steps that can cover approximately 
eight orders of magnitude. The time steps are determined 
by the local Courant-Friedrichs-Lewy criterion l|Monaghanl 
1989) ensuring that information is not carried over length 
scales larger than the local kernel size that is defined by the 
20 nearest neighbor particles in two dimensions. In addition 
we require that the fractional change of the kernel size of 
each particle is limited to 10 percent. Thereby we ensure 
that particles require several time steps to pass through the 
interface of a strong shock. 

For the fixed-grid simula tions we used th e gas- 



kinetic grid code Proteus p^rendergast fc Xiil 19931 : 
ISlvz fc Prendergast|[l999l : IXull200ll : IHeitsch et al.ll2004h in 
an ide ntical implementation as described by Heitsch et al.l 
(|2008l ). Proteus allows the explicit control of vis cosity and 
heat conductivity at fixed spatial scales (see ISlvz et al] 
|2002| . for a detailed description). 



2.1 Heating and Cooling 

To compute the heating and cooling rates we have used the 
rates for optically thin atomic lines from W olfire et al.l (|l995l , 
l2003h . The cooling curve covers a range of 10"^ sS n s; 
10^ cm"^ in density and 30K s: T sC 2 5 x 10"^ K in temper - 
ature and is identical to that used by IHeitsch et al.l l|2006l ). 
Dust extinction becomes important above column densities 
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of A''(HI) ^ 1.2 X lO^^cm"^, which are reached only in the 
densest regions modeled. Thus, we use the unattenuated UV 
radiation field for grain heating, expecting substantial un- 
certainties in cooling rates only for the densest regions. The 
ionization degree is derived from a balance between ioniza- 
tion by cosmic rays and recombinatio n, assuming that Ly a 
photons are directly reabsorbed. As in lHeitsch et al.l l|2006h . 
heating and cooling is implemented iteratively as a source 
term for the internal energy e of the form 



J = nr(r) - n^A(r) 



(1) 



in units of energy per volume per time. Here, F is the heat- 
ing contribution (mainly photo-electric heating from grains) , 
and nA is the cooling contribution (mainly due to the ClI 
line at 158/im). To speed up the calculations, equation ((T| 
is tabulated on a 2048^ grid in density and temperature. 
For each cell and iteration, the actual energy change is then 
bilinearly interpolated from this grid. 

The implementation of heating and cooling in VINE 
follows closely that for Proteus described above. For each 
particle i the change of the specific internal energy, u, in 
units of energy per unit mass per unit time is computed as 



=n.F(rO-n?A(r,) 



(2) 



If the current integration time-step of the particle is larger 
than one percent of the local cooling time the cooling is 
computed iteratively in steps of one percent of the local 
cooling time. Using periodic boxes of constant density and 
temperature at different values we have tested that this im- 
plementation results in the correct cooling rates. 



2.2 Initial and Boundary conditions 

Two opposing, uniform, identical flows in the x-y compu- 
tational plane initially collide head-on at a sinusoidal inter- 
fa ce with wave numbe r ky = 1 and amplitude A (see Fig. 1 
m iHeitsch et allbood '). The incoming flows are in thermal 
equilibrium at a density of no — 1.0 cm~^, a temperature 
of To = 8.5 X 10^ K and a velocity of 21 km s~^, typical for 
parameters of the warm neutral medium. The box length is 
44 pc. 

We ran identical model sets with VINE (model names 
headed by "V") and Proteus (names headed by "P"). The 
model name V0256 specifies an SPH model with initially 
256'^ particles, while e.g. P2048 denotes a grid model with 
2048^ cells. 

Initially we distribute the SPH particles on a grid in 
the square domain. The particle mass resolution is 7.2 x 
10~* M© pc~^. We create new particle layers at the x- 
boundaries at the initial grid resolution with the appropriate 
velocities. Thus, the total particle number of the simulation 
is continuously increasing, starting with with 256^ — 65536 
particles. At 9 Myr, there are already 640, 000 particles in 
model V0256. This leads to a significant slowing down as 
the simulation proceeds. 

We use p e riodic boundary conditions in the y-direction. 
iHeitsch et all (|2006h explored the effect of the choice of 
boundary conditions on the mass budget and the morphol- 
ogy of the forming clouds. They found that although open 
boundaries allow a small fraction of the gas to leave the 



simulation domain (or active region), the total mass bud- 
get and the morphology is only slightly affected. The cool- 
ing timescale is much shorter than the flow timescale and 
the sound crossing time, so that the gas is compressed and 
cooled down before it can feel the boundaries. Only gas close 
to the boundaries will be affected. 



3 RESULTS 

3.1 Morphologies 

We begin with the flow morphologies (Fig. [T]), comparing 
density and temperature maps of models V0256 (left) and 
P1024 (right). Densities and temperatures have been scaled 
to the same ranges, and for both codes, the central quar- 
ter (22^ pc^) of the simulation domain is shown at times 
t = 2.3, 4.2 and 9.2 Myr after initial flow contact. The most 
noticeable difference is the much more filigree structure of 
the cold gas exhibited in the SPH simulations. Comparing 
the initial resolution (256^ vs 1024^) against the effective 
resolution indicates the resolving power of the Lagrangian 
method. The smallest smoothing lengths in model V0256 
correspond to 2 x 10""^ pc, while the grid resolution of model 
P1024 is 4 X 10"^ pc, and that of P8192 5.3 x 10"^ pc. 

The SPH models tend to show a broader post-shock 
region (greenish colors) than the grid models. The densities 
for V0256 and P1024 are scaled identically (between —0.7 < 
logn < 3.5), as well as the temperatures (between 1.0 < 
logT < 4.4). V0256 clearly shows thinner filaments, but at 
approximately the same density range, as model P1024 - a 
consequence of the density-dependent spatial resolution of 
SPH. 

Figure [2] shows a zoom sequence of models V0256 (left) 
and P8192 (right). The spatial extent of each zoom level is 
indicated in the lower left of the respective panel. The spa- 
tial resolution of model P8192 (5.3 x 10~^ pc) is only a factor 
of 2 away from resolving the Field leng th ([F ield 196 f^) for a 
realistic value of the heat conductivity (|Kovama fc Inutsukal 
[2004). A closer study of the corresponding temperature 
maps (not shown) indicates that the cooling length espe- 
cially between the warm ambient and cold dense gas is re- 
solved in P8192: the cold dense regions start to show "halos" 
of intermediate temperature material. 

Comparing V0256 and P8192 at the largest scale (L = 
19 pc), it is already obvious that the SPH model reaches 
higher densities, at a more filamentary appearance. Con- 
versely, stuctures in P8192 seem to be more extended, and 
less filamentary. This changes when moving to smaller scales 
(down to 2 pc). V0256 exhibits mostly coherent filaments at 
high densities (n > 10^ cm~"^), while P8192 shows a whole 
population of small cloudlets including turbulent wakes. 
These are conspicuously absent in V0256. Thus, while the 
density-dependent spatial resolution allows for higher den- 
sity contrasts and smaller structures in V0256, the constant 
spatial resolution in P8192 leads to a better tracking of small 
low-density structures. The differences suggest (compared to 
the density ranges in Fig. [T]) that VINE reaches an effective 
grid resolution in the dense gas of 1024^ for (initially) 256^^ 
particles. The highest zoom also suggests that the density 
contrasts are more sharply resolved in P8192 than in V0256, 
leading to smoother gradients in the SPH model. 
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Figure 1. Morphologies of colliding flows. Shown arc logarithmic density and temperature maps of the central quarter of the simulation 
domain, at times t = 2.3,4.2 and 9.2 Myr, left for model V0256, and right for model P1024. The effective spatial resolution of V0256 
leads to more filigree structure than in model P1024, with fixed resolution. Densities are scaled identically between —0.7 < logn < 3.5 
for both models, and temperatures between 1.0 < logT < 4.4. 



The highly structured "slab" of model P8192 and a 
comparison with P1024 shown in Figure [T] emphasizes the 
importance of resolution for the problem of turbulent ther- 
mal fragmentation: only at high Reynolds numbers, a tur- 
bulent cascade can develop. Proteus runs with a specified 
viscosity, measuring u = 8.9 x 10^^ cm^ s~^ for our parame- 
ters. Using half the period of the geometric perturbation in 
the collision interface as typical length scale, and the inflow 
velocity as typical flow velocity, results in a formal Reynolds 
number of Re ^ 1.6 x 10^ for P1024, and Re « 1.3 x 10'' 
for P8192. Since the condensation mode of the TI responsi- 
ble for the fragmenta tion grows first on the smallest scales 
l|Burkert fc Linll2000l ). the generation of small-scale struc- 
tures due to turbulence on which the TI can feed will be 
crucial for the overall evolution of the system. 

3.2 Mass Fractions 

Figure [3] shows the mass per length for the cold thermally 
stable regime (T < 300 K), the warm thermally unstable 
(300 < T < 3000 K), and the warm thermally stable regime 
(3000 K< T), as a function of time. For both codes, the mass 



of cold gas increases nearly linearly with time once the ther- 
mal instability has set in. Thus, despite all its substructure, 
the interaction zone acts like a one-dimensional slab in terms 
of mass collection. The slight offset in time between the SPH 
and grid models is caused by the stronger effect of the ar- 
tificial (bulk) viscosity in SPH to prevent the penetration 
of particle layers. The gas thus reaches higher temperatures 
and requires more time to arrive at densities beneficial for 
strong cooling (and mass accumulation). This leads to a con- 
stant mass offset of ~ 70 Mq pc~^ in the cold gas between 
Proteus and VINE models, corresponding to ~ 20 Mq pc~^ 
less intermediate temperature gas and ~ 50 Mq pc~^ less 
high temperature gas in the Proteus models. Mass histories 
of models P4096 and P8192 are indistinguishable. 

While the SPH models show less mass in the cold tem- 
perature regime, there is more mass in the intermediate and 
warm regime. Both codes evolve the warm and intermedi- 
ate regime qualitatively similarly: the mass in both regimes 
stays approximately constant, since they are only transitory 
stages of the gas on its way to the cold regime. 

Considering the widely differing numerical methods, the 
mass fractions agree surprisingly well. Especially the prob- 
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Figure 3. Total mass per length in the three temperature regimes T < 300 K, 300 < T < 3000 K and 3000 K < T as a function of time. 
Solid thick lines denote the grid models P2048, P4096 and P8192, while dashed line shows model V0256. Mass histories for models P4096 
and P8192 are indistinguishable. The constand mass-offset between VINE and Proteus models is due to the higher artificial viscosity in 
VINE. 



lem of over-cooling in SPH-methods l|Pearce et al.l \l99^ ) 
does not seem to dominate the mass budget. 

3.3 Clump Mass Function 

The strong fragmentation due to the thermal instability 
begs the question what the resulting clump mass function 
(ClMPfl of the cold clouds will be. While the detailed phys- 
ical processes linking the CIMF to the core mass function, 



ter of debate (e.s. Klessen & Burkert 20001: Williams et al. 


200d 


iPadoan & Nordlund 20021: Ballesteros-Paredes et al. 


2000 


lElmes:reen"2007': iHennebelle & Audit 12007: Dib et al. 


2008 


Smith ct al. 2009|), a model of rapid star formation 



in colliding flows should explain what physical processes 
are setting the CIMF during the formation of the molec- 
ular cloud. The two main fragmentation agents are cool- 
ing and self-gravity, the latter of which we have not imple- 
mented in the current models. While we envisage self-gravity 
to be important on global sca les (|Burkert fc HartmannI 
l2004l : iHartmann fc Burkertll2007l) for finite clouds (i.e. non- 
periodic boundaries), on local scales, the gas has to cool 
substantially before the Jeans lengths get small enough for 
gravity to dominate. Which mechanism then dominates set- 
ting the CIMF? Because of the early and strong fragmen- 
tation due to thermal (and turbulent) effects, the CIMF 
mig ht well be determined e arly on in the cloud evolution (see 
also I Hennebelle et al.ll2007l for 2D models, and lHeitsch et ahl 
l2008l for 3D models'). 

The main strength of VINE over Proteus is of course 
the adaptive spatial resolution. To identify clumps of cold 
gas in model V0256 we have used a friends-of-friends algo- 
rithm with a linking length of 0.04 pc, which is about ten 
times the minimum kernel softening length at this time, and 
a minimum of 20 particles per clump. In total we found 418 



^ We use the abbreviation CIMF to distinguish between clump 
and core mass function, the latter of which is usually abbreviated 
as CMF in the literature. 



individual clumps of gas with temperatures T < 300 K. The 
CIMF is plotted in Fig. H at a model time of 9.2 Myr. For 
comparison we have computed the CIMF of model P2048, 
P4096 and P8192 using the same method. Instead of parti- 
cle positions, we used the cell positions of the cold gas as 
input for the friends-of-friends algorithm, again with a link- 
ing length of 0.04 pc. For the Proteus models, the linking 
length is a factor 2 to 8 larger than the (linear) cell size. 
As expected, the number of identified clumps increases for 
the grid simulation with resolution (153 clumps for P2048, 
495 for P4096, and 860 for P8192). At high masses, P8192 
and V0256 agree, whereas there is a clear deficit in clumps 
at lower masses for the grid models, with a factor of three 
difference between V0256 and P8192 at M 1 Mq pc"\ 

A direct comparison of our CIMF to observations is 
problematic, for two reasons. First, the restriction to two 
dimensions will emphasize compression (cooling) over shear 
flows and vorticity (turbulence) since the gas flows are 
more compressible than in three dimensions (see Fig. 13 of 
iHeitsch et al.ll200^ . Second - as the referee pointed out -, 
we chose to calculate ClMFs at a time (9.2 Myr) at which 
a sufficient amount of gas will have accumulated for self- 
gravity to affect the evolution of th e cloud. This can be 
seen c omparing models Hfl and Gfl of lHeitsch &: HartmannI 
(|2008l . see their Fig. 2) at (their) 9.9 Myr. We chose this 
rather late time for our two-dimensional models, since the 
suppression of the third dimension leads to rather coherent 
sheets filaments still imprinted by the initial conditions at 
earlier times. 

Keeping the above caveats in mind, the SPH slope is 
consi stent with the clu mp mass slopes fou nd by (|Li et al.l 
l2007l , see also Fig. 23 of iRagan et"al]|2009ll for Orion, and 
with earlier me a surem ents by Heithausen et all l|l998l l or 
ISchneider et al.l (|2002l l. Overall, the slopes of both meth- 
ods are shallower than those observed for low-mass star 
forming regions such as Perseus l|Curtis fc Richeij|2010h . Yet 
these authors point out that such differences might well be 
caused by the way in which clumps are defined, since dif- 
ferent clump-finding algorithms can result in very different 
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Figure 2. Zoom sequence of V0256 (loft) and P8192 (right), log- 
arithm of density, shown at t = 9.2 Myr, at (from top to bottom) 
a scale of 19, 9, 4.25, and 2 pc. Color scales are identical for both 
models, highlighting the differences. The full simulation domain 
would measure 44 pc a side. The smallest dense structures (at 
a zoom of 2.0 pc) have a size of < 0.1 pc. Turbulent wakes are 
clearly visible in model P8192. 
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Figure 4. Histogram of clump masses (CIMF) for VINE 
(V0256,dashed Une) and Proteus (P2048, P4096, P8192, solid 
lines). Sit t = 9.2 Myr. The dotted line shows a slope of —0.85 
as determined for massive clumps in the Orion Molecular Cloud 
l|Li et al.ll2007l 'l. 
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slopes. In any case, for the reasons given above, our ar- 
gument here focuses more on the comparison between the 
numerical methods (for which we use the identical clump- 
finding mechanism), rather than on a detailed comparison 
with observations. 

3.4 Line-of-sight Velocity Distributions 

In Fig [5] we show the mass-weighted line-of-sight velocity 
distribution of the cold gas (T < 300 K) taken along the 



Figure 5. Linc-of-sight velocity distribution of the cold gas (T < 
300K) along the inflow direction. Dashed line: VINE, solid line; 
Proteus. 



inflow direction. The velocity dispersion is « 6.8 km s~^ for 
VINE and f» 8.5 km s~^ for Proteus. The sound speed of 
the cold gas is ~ 0.75 km s^^ at T = 40 K, i.e. the cold 
gas motions are nominally supersonic. 

While the overall agreement of the distributions is 
reassuring, the differences are caused by resolution ef- 
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Figure 6. Ratio of internal velocity dispersion and internal sound 
speed (i.e. the internal Mach number) for all cold clumps. For all 
models, most of the clumps have subsonic internal motions. 



Figure 7. Internal sound speed of the cold clumps (T < 300K). 
The clumps found in model V0256 have systematically lower in- 
ternal dispersions with a smaller spread. 



fects: At higher reso lution, smaller scales start to thermally 
fragment first (e.g. iBurkert fc Liiil I2OO0I ) suppressing the 
growth of the dynamical instabilities (e.g. IVietri et al.lll997l : 
[Hucckstacdt 2003). Since the SPH models have higher spa- 
tial resolution in their evolved state, we expect them to 
be more dominated by thermal effects rather than by the 
dynamical instabilities, thus leading to a (slight ly) smaller 
velocity dispersion (see also iHeitsch et aLll2005h . Also, we 
expect the generally higher dissipative nature of S PIf com- 
pared to grid methods (e.g. iKitsionas et al.ll2009l ) to con- 
tribute. 

With the sound speed in the cold gas Cs ~ 0.75 km s~^, 
the internal velocity dispersion of the cold clumps turn out 
to be mostly subsonic (Fig. This confirms earlier find- 
ings from numerical simulations (iKovama fc Inutsukall2002l : 
I Audit fc HennebelldliooBI : iHeitsch et al.l l200d). aiidshould 
not come as a surprise in view of the formation mecha- 
nism of the cold dense clumps. The Proteus models have 
a broader distribution of internal velocity dispersions than 
the VINE models, and the peaks shift to smaller veloci- 
ties with increasing resolution, indicating that the thermal 
length scales ar e not completely resolved at the lowest reso- 
lutions (see also lHennebelle fc Auditl2007l ). The distribution 
for the V0256 clumps is much narrower and peaks around 
0.25 km s"\ For comparison we show the distribution of 
the internal sound speeds of the individual clumps in Fig. [T] 
Both distributions peak around 0.45 - 0.50 km s^^. However, 
the clumps in the SPH simulations tend to cluster at lower 
temperature (T < 100 K) and there is less gas at higher 
temperatures than in the grid simulations (Fig. [S]). This is 
probably a result of SPH intrinsic over-cooling. 



0.30 _ 




Temperature 



Figure 8. Fractional mass as a function of temperature in the 
cold phase (not limited to clumps). There is relatively more low 
temperature (T < 40 K) gas in model V0256 than in the Proteus 
models. 

4 SUMMARY AND CONCLUSIONS 

Lagrangian hydrodynamical methods such as Smoothed 
Particle Hydrodynamics (SPH) are an attractive tool to sim- 
ulate astrophysical problems with a wide range of spatial 
scales, such as self-gravitating objects or thermally unstable 
flows. However, due to the numerical implementations, the 
reliability of SPH for such applications has been questioned. 
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We compare a SPH implem entation (VINE, I Nelson et alj 
I2OO9I : Wetzstein et al.ll2009l'l and a grid-based method (Pro- 
teus , Prcndcrgast fc Xu"l993l: ISlvz fc Prendergast|[l999l : IXul 
l2001i : Heit sch et al., 2004 ) in an application to flow-driven 
molecular cloud formation. This is a challenging astrophys- 
ical problem for any numerical method, because of the tur- 
bulence, the high density and temperature contrasts, and 
the rapidly shrinking scales of the cold clouds. 

The morphologies between particle and grid based 
methods agree well globally, yet the cold gas structures are 
more fragmented in the SPH models (Fig. Temperature 
and density gradients in SPH tend to be "rounder" than in 
the grid method as a consequence of deriving physical quan- 
tities by averaging over a set of nearest neighbors. However, 
this drawback is countered by the higher spatial resolution, 
yielding more filigree structure in the SPH models. 

The mass fractions of gas in different temperature 
regimes differ slightly (Fig. [3]). SPH models start collecting 
mass in the cold stable regime later than the grid models, 
leading to a constant offset in the mass with time. This 
is a consequence of the strong overheating in the initial 
flow collision due to SPH's artificial viscosity. The slopes 
of the ClMFs (Fig. |4| are consistent with res ults from ear- 
lier numerical models (|Hennebelle et aLll2007l ) with similar 
physics. They are also consistent with observational re sults 
i|Schneider et al.ll2002l : lLi et al.ll2007l : lRagan et al]|2009l ). yet 
such a comparison should not be over-interpreted because 
of the restricted dimensionality of our models as well as 
the fact that self-gravity will have affected the evolution 
of the cloud at the time w e chose to analyze the models 
ijHeitsch fc Hartmannll2008h . 

The line-of-sight velocity distributions have similar 
widths (Fig. [5}, namely « 6.8 km s"^ for VINE and 
~ 8.5 km s~^ for Proteus, rendering the gas nominally 
supersonic, while both methods find that the internal 
velo city dispersion of the cold clumps is subsonic (see 



also iKovama fc Inutsukal l2002l : lAudit fc Hennebelld I2OO5I : 
iHeitsch et al.ll2006l l. The slightly narrower distribution for 
SPH is a consequence of the higher resolution, emphasizing 
the thermal instabilities over the dynamical ones. 



The effects of overcooling ijPearce et al.l [l999h do not 
seem to affect the overall mass budget (Fig. [3| within the 
cold stable regime, but they do manifest themselves in the 
low-temperature tail of that phase (Figs. [7] and [S]). This sug- 
gests that the TI controlling the transition from the warm 
stable to the cold stable regime is not strongly affected. 

Our results suggest that the SPH implementation in 
VINE is suitable to model thermally bistable, dynamical en- 
vironments such as occurring in the flow-driven formation of 
molecular clouds. Yet, quantitative predictions will slightly 
differ between particle and grid methods, with a bias of the 
VINE results to lower temperatures and smaller fragments 
in the cold gas, but a lower total mass in the cold, thermally 
stable regime compared to Proteus. In terms of mass spec- 
tra, an initial particle resolution of 256^ is (at later stages) 
equivalent or better than a fixed grid of 4096^ cells, with 
promising consequences for the determination of e.g. clump 
mass functions. 
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